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The two-dimensional (2D) quantum Heisenberg antiferromagnet at zero 
temperature, with locahy frustrating magnetic defects, is studied. We con- 
sider two types of defects - an isolated ferromagnetic bond (FMB) and a 
quantum impurity spin, coupled symmetrically to the two sub-lattices. A 
technique is developed to study strong defect-environment interaction. In the 
case of a FMB we find that, contrary to the previous linear spin- wave result, 
the local magnetization stays finite even for strong frustration. For an anti- 
ferromagnetic coupling of a quantum impurity spin with its environment, we 
find a local change in the ground state. All our calculations are compared 

with numerical results from exact diagonalization. 
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I. INTRODUCTION 



The two-dimensional quantum Heisenberg antiferromagnet (HAFM) has attracted a lot 
of attention in the last several years, mainly because of its relevance to the physics of the 
high-Tc material^. At zero hole doping, these compounds form planes of localized Cu spins 
and their physics is well described by the spin-1/2 square lattice HAFM. There is strong 
numerical evidence that this system has long-range order (LRO) in its ground state at 
T = 0. Due to quantum fluctuations, the staggered magnetic moment was found to be 
reduced to about 61% of the classical valueB. Doping with excess oxygen, which is crucial 
for superconductivity, creates spin-1/2 holes on the oxygen sites, located in between the Cu 
sites. The holes interact magnetically with their Cu neighbors. This interaction leads to 
frustration of the Cu spins and ultimately to destruction of magnetic LRO, above a critical 
dopant concentration. 

One way to describe the doping process is via an effective one-band model, like, e.g. the 
t — J mo del, in which the holes hop in the antiferromagnetic (AFM) environment. 

However, the extreme limit of static holes, which frustrate the AFM order only locally, 
is also believed to have relevance. The basic idea, put forward by Aharony et a/.i, is that a 
hole, interacting with its two neighbors creates an effective ferromagnetic interaction between 
them. The magnetic properties then are described by a HAFM with a given concentration 
of such ferromagnetic bonds (FMB). As a start, it is important to understand the effect 
that one isolated FMB has on the AFM environment. This problem was studied by Lee and 
Schlottmanni and Aristov and Maleevi in the linear spin- wave approximation (LSWA). They 
found that the magnetization of the two spins, connected by the bond decreases dramatically 
as the strength of the bond increases, and noted an instability that occurs in the theory at 
the "classical" transition point (the point where a local triplet formation is favorable without 
taking quantum fluctuation into account). There have been several extensions of this work, 
taking into account finite temperature and finite concentration of bondsii. One can also 
treat the frustration effects more realistically, by considering a model of frustration, in which 
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a static hole interacts with its two neighbor spins. Oitmaa, Betts, and AydinO studied this 
model numerically, via an exact diagonalization on small clusters, and found a discontinuous 
local phase transition for an antiferromagnetic hole-neighbor interaction, whereas no change 
in the ground state occurred for a ferromagnetic coupling. A similar model was studied by 
Clarke, Giamarchi, and Shraimani in one dimension via bosonization. 

In the present work we study the two models with local frustration, mentioned above - an 
isolated FMB in the 2D HAFM, and a quantum defect spin, coupled symmetrically to the 
two sub-lattices. Our motivation is twofold. First, the basic works on the FMB problemBi 
use the LSWA, which is only valid for small FMB strength. In the original frustration model 
however, Aharony et a/.i concluded that the induced FMB coupling should be large. Thus, it 
is important to study the model for strong coupling, and, along the way, assess the region of 
validity of LSWA. Second, the two-dimensional spin defect model, which provides a realistic 
frustration mechanism, has not been studied analytically, as far as we know. The present 
work, thus, extends the previous work of Oitmaa, Betts, and Aydini on this problem. 

Our analytical approach is specifically designed to treat the strong coupling regime. The 
main idea is that since the perturbations are local, one should first solve the defect problem 
separately, and then take into account the defect-AFM environment interaction. For the 
FMB problem this amounts to considering the two-body Green's function of the two spins, 
connected by the bond. Analogously, for the spin defect problem, we introduce a three-body 
Green's function (of the spin defect and its two neighbors). All our analytical results are 
compared with numerical simulations, based on exact diagonalization on small clusters. A 
short account of the main results of this work appeared previously in Ref. [10]. The rest 
of the paper is organized as follows. Sec. II describes the basic idea behind our theoretical 
approach and, as a first application, we consider the HAFM with no defects. In Sec. Ill 
we present our analytical, as well as numerical solution of the FMB problem. Sec. IV deals 
with the spin defect problem from the same perspective. Our conclusions are summarized 
in Sec.V. 



3 



II. METHOD DESCRIPTION. 



Consider the two-dimensional spin-1/2 Heisenberg model on a square lattice at zero tem- 
perature: 

Hj = JY. Sj,. (1) 

<hj> 

The summation is over nearest neighbors only, and J > 0. We set J = 1 from now on. It 
is knownB that there is LRO in the ground state of Hj and the staggered magnetic moment 
m =< Sf >= 0.303. 

Now assume we introduce additional interactions (of typical strength K) between two 
nearest neighbor spins. Let us call such a configuration a spin defect. We have in mind, for 
example, a FMB, connecting two of the spins. Since such a defect is local, it does not affect 
significantly the LRO far away from it. Therefore, in order to treat the strong coupling 
regime {K ~ J), one must first take into account the interaction between the spins in the 
defect and then treat the defect-AFM environment interaction. 

As a simple illustration of the above idea consider the HAFM (1) without any additional 
interactions. Let a "defect" consist of only one spin 5*1, and separate (1) into an AFM 
background (spin-wave) Hamiltonian and a defect-neighbor interaction: 



H = H,^ + S^,j2S! + llstj:S-+ h.c.j, 

1=2 I ^ i=2 J 



(2) 



i^s«, = E^k(4«k + 4/5k). (3) 

k 

Here is /?k and ctk are the usual spin-wave operators with dispersion0: 



= 2y 1 - 7^, 7k = -(cos(A;^.) + cos{ky)). (4) 

Next, the Ising part of the interaction in (2) is taken into account in the mean- field 
approximation, i.e. we replace: J2^=2 4m. For the spin 5*1 it is convenient to use the 

fermion representation: 
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st=^^^, si=-m^^~¥^^^), (5) 

where and create (acting on their common vacuum) a component 1/2 and —1/2 
respectively. In order for the two fermions to represent a spin-1/2 operator, they have to 
satisfy the constraint: + ^j^E'; = 1, i.e. the physical states have only one fermion. In 

this representation the Hamiltonian becomes: 

H = Hsw + Ho + Hint, (6) 



i/o = 2m(vl>Ivl/^-vl/|vl/^), (7) 

/ g 

Hint = \ T7^{^i E 7k(Mk/?k + t^kftlk) + h.c, (8) 
V k 

where N is the total number of lattice sites, and = '^^ = — sign(7k)y^— | + ^ 

are the parameters of the Bogoliubov transformation. Here we have assumed that the spins 
Si,i = 2, .., 5 belong to the B sub-lattice (spin down). 

We define the Green's functions for the two fermion and spin-wave species: 
G^{t) = -i< r(^^(t)^i(0)) >, ^ =T, i, D%t) = -i< T{Pk{t)Pl{0)) > (and similarly 
for ttk)- The unperturbed Green's functions {Hint = 0) are: 

Gui^)= , D{k,u;) = (9) 

uj ± 2m + 10 uj — e\^ + Id 

Next, the term Hint is treated in perturbation theory. The self-energy to one-loop order is 
given by the diagram in Fig. la : 

Thus the spectral function of the particles to this order iJl3: 

A^{e) = --ImG:"*(£) = 6{e + 2m - S|(e)) = Z6{e - e*), (11) 

TT 

where e* is the solution of the equation: 



£* + 2m-ET(£*) = 0. (12) 
The renormalization factor Z is determined from (11): 

The ground state wave function of the ^ particle - spin-wave system is then: 

10= V^*||0>+^Ak*;4|0>, (14) 

k 

where |0 > is the vacuum for the and the spin- wave operators: 

*^|0 >= ak|0 >= /?k|0 >= 0. (15) 

From now on we will also adopt the short-hand notation: 

|/x>=*i|0>, |/x,/5>=^'J,4|0>, >= *i4|0 >, ix=\A- (16) 

The wave function (14) is a sum of coherent (first term) and incoherent part. In the coherent 
part, \fZ represents the probability amplitude of having the particle ^| (i.e. the spin 5*1 
being up) in the ground state of H. The incoherent part comes from the admixture of 
intermediate states, each of them entering with probability amplitude 

A = ^^-^l^^-'l (17) 
e — 2m — £k 

It is easy to see that: 

$:iAkr=i-z, (18) 

k 

which reflects the normalization condition for \G >. Finally, the magnetization is: 

<G\Sl\G>^)^Z-h^\-Z)^^.T(. (19) 

The agreement between (19) and the spin-wave theory result (m = 0.3) is quite good. 
There are two ways to improve our calculation. First, two-loop (see Fig. 2b), and higher 
order corrections can be included. Let us note that due to spin conservation there are no 
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vertex corrections to the self-energy (of the type, shown in Fig. 2c). A similar phenomenon 
occurs in the t — J mo deli. Second, we could define clusters of more than one spin as our 
"defect". Although the calculation accuracy is expected to improve by taking the above 
two points into account, the leading order result is quite reliable, and, for the problems 
considered in the next two sections, is in good agreement with numerical simulations. 



III. FERROMAGNETIC BOND PROBLEM. 
A. Theory. 

Consider a ferromagnetic bond of strength K, connecting the sites 1 and 2: 

H = Hj- KS1.S2, (20) 

where i^' > 0. The second term in (20) frustrates the LRO in the neighborhood of the defect 
and therefore leads to local suppression of the magnetization. 

Following the idea, outlined in Sec. II, we separate the interactions into a part within the 
defect (i.e. between the spins Si and 5*2) and a defect - spin-wave interaction. The Ising 
part of the latter is treated in mean-field approximation, whereas the transverse part {Hint 
below) is our perturbation: 

H = Hsw + Ho + Hint-, (21) 
HQ = -{K-l)Si.S2 + ?>m{S'^-Sl), (22) 

H^nt = E r(k)e*'=^(Mk4 + ^k/3-k) + ^1+ E r(k)(Mk/3k + t^k«Lk) + h.c.j , (23) 

where 

r(k) = 2cos(A;j^) + exp {ik^), (24) 

and Hsw is given by (3). In (22) and (23) it is assumed that Si belongs to sub-lattice A 
(spin up) and 5*2 - to sub-lattice B. 



We introduce a fermion representation for the spins Si and S2, similarly to (5), and call 
the corresponding fermions ^-f^^ and ^1,1, respectively. Thus (22) and (23) should be thought 
of as expressed in terms of these fermions. In order to diagonalize Hq it is convenient to 
define the two-particle matrix Green's function Gnu' 



^ < r(*T(t)$^(t)vi/|$|) > < T(*T(t)$^(t)vi/|${) > ^ 



\^ < r(*i(i)$tW*{$I) > < r(*i(i)$TW*I*T) > 

As (25) suggests, the diagonal elements (11 and 22) correspond to the two-particle states 
|1 >= I tJ,> and 1 2 >= | J,t>, respectively, where the first arrow represents the spin 5*1 and 

— * 

the second one - the spin 82- The off-diagonal components represent transitions between 
these two states. For future purposes it is convenient to define also the states: |3 >= | tT> 
and |4 >= | J,J,>. The unperturbed Green's functions (corresponding to Hq) are: 

g....M- ^_(;,_lV4±3m + « - G..M-G..M- ,^(^_\)^,^,, . (26) 
Next, we evaluate the self-energies to lowest order in perturbation theory with respect to 
Hint- It is, in fact, more convenient to use the Releigh-Schroedinger perturbation theory, 
rather than Feynman's diagrammatic technique. Using notations, similar to (16), the ex- 
pression for Ell is: 

^^^^^)-^.+ (X-l)/4-., + .+ (X-l)/4-.,- ^'^^ 
Here by Hint we mean the interaction Hamiltonian at fixed wave vector, i.e. (23) without 
the k summation. Evaluating (27), we get: 

^''^'^-N\s + iK-l)/A-s^- ^^^^ 



(25) 



Similar calculations for the other two Green's functions give: 

1 ^ \T{k)\^vl ^ , , 1 Re[r2(k)e 



This is the one-loop result which is sufficient for our purposes. Let us mention that higher 
loop orders of perturbation theory for G^i, are more complicated than those for the one- 
particle propagator of Sec.II, because vertex corrections are allowed. 



0. (30) 



In order to evaluate various correlation functions, we proceed similarly to the calculation 
of the previous section. Since the Hamiltonian mixes the fermionic states |1 > and |2 >, the 
equation for the effective energy level £*, corresponding to (12) now has the form: 

e*-{K- l)/4 + 3m - ^u{e*) -{K - l)/2 + Ei2(£*) 

-{K - l)/2 + ^^2is*) e*-{K- l)/4 - 3m - H^^ie*) 
The correct normalized eigenstate is of the form |12 >= > +/X2|2 >, where {/-ti, H2) is 
an eigenvector of the matrix in (30), and /if + /i^ — 1. 
The ground state wave function can be written as: 

IG >= v^l 12 > + incoherent part, (31) 

where the incoherent part represents the contribution of the intermediate states |3 > and 
|4 >. The normahzation factor Z is defined as: 

where the self energy E: 

t{e) = nl^uie) + //^E22(£) + 2niH2^uie). (33) 
The average of any spin operator O in the state |G > is: 

< > = < .101. > z . _ s <'\o\i> 

^< 12\d\12 > Z+^^^-^ J2 <i\d\i> . (34) 

2 j=3,4 

Using the above formula, we get for the magnetization and the longitudinal and transverse 
spin-spin correlation functions: 

M =< St >= d^z, (35) 
C^l, 2) ^< SlSl >= -\ - \{Z - 1), (36) 

Ct{\ 2) = ^ < 5+^2" + h.c. >= ii^ii^Z. (37) 

The numerical evaluation of the sums (28) and (29) and the solution of (30) are straightfor- 
ward. The correlators are then computed by using Eqs. (32,33,35-37). 



B. Exact diagonalization studies. 



To estimate the accuracy of our analytic approach we have obtained numerical results by 
exact diagonalization of small clusters. Instead of using periodic boundary conditions, we 
have applied a staggered magnetic field in the z direction to spins on the boundary of the 
cluster. This breaks the sublattice symmetry and allows us to compute single-spin averages 
as well as to distinguish between longitudinal and transverse correlations. 

This procedure of course also breaks the translational invariance of the cluster and will, 
in general, lead to increased finite size effects and slow convergence to the bulk limit. We 
have chosen a cluster of N = 18 sites and a boundary field to give < >= 0.3 on the 
boundary spins. Extremely good qualitative agreement was obtained with the analytical 
results, and tuning of the field would achieve even better agreement. This however is not 
warranted without, at the same time, a systematic extrapolation of the cluster results to 
the thermodynamic limit. 



C. Results. 

The results obtained by the Green's function method as well as the exact diagonalization 
results are summarized in Figures 2 and 3. 

Figure 2 plots the local magnetization M and the total spin-spin correlation function 
across the FMB C(l,2) =< S1.S2 >= Ci(l,2) + Ct(1,2) as a function of frustration K. 
We have also calculated M by using LSWA, essentially following Lee and Schlottmanni. 
By treating as a perturbation one can, in fact, exactly solve the Dyson equation for 
the spin-wave Green's functioni. The LSWA result shows that the magnetization vanishes 
around K = 1.9 and an instability occurs at = 2, which is the point when a local triplet 
formation is expected to occur classically. 

Our result (35) shows qualitatively different behavior. Up to K = 1 the two curves 
follow each other closely, both predicting a slight increase of M at that value (corresponding 
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to a missing bond). A similar suppression of quantum fluctuations was observed in the 
case of a missing siteEl. However, beyond = 1 we predict that the local magnetization 
decreases slowly and stays non-zero even for large K. On the other hand, the correlator 
C(l,2) changes sign and becomes ferromagnetic at around K = 2.1. As shown in Figure 3, 
it is, in fact, the transverse part of the correlator that changes sign, while the Ising part has 
antiferromagnetic sign for all K. With increasing K the transverse part increases and the 
Ising one decreases (in magnitude), but remains non-zero, which is consistent with a finite 
magnetization. 

The exact diagonalization results agree qualitatively well with those obtained by the 
Green's function method. 

We would like to stress that, as our results show, the range of validity of the LSWA is 
limited to small frustration only {K < 1). Beyond this point the number of generated spin 
waves is large and the interactions between them become important. The vanishing of the 
magnetization in LSWA is thus an artifact of the approximation. Our treatment, on the 
other hand, is applicable far beyond K = 1. 



IV. IMPURITY SPIN PROBLEM. 

In this section we consider local frustration due to an additional quantum spin 
a (cr = 1/2), coupled symmetrically to the two sub-lattices: 

H = Hj + La.{Si + S2)-Si.S2. (38) 

For either sign of the coupling L, the second term in (38) frustrates the spins 5*1 and 5*2. It is 
also assumed that the impurity spin effectively removes the superexchange between themil. 

Similarly to the approach used in Sec. HI, first we diagonalize the second term in (38) 
exactly, and then treat the interaction of Si and 5*2 with the AFM background (spin waves) 
perturbatively. Thus we have: 

H = Hgyj + Hq + Hint, (39) 
11 



where 

Ho = La.{S^ + ,^2) + 3m(5| - SI) , (40) 

and Hgw and H^nt are given by (3) and (23) respectively. We have assumed, as in the 
derivation of (23), that Si e (sub-lattice A), S2 & (sub-lattice B). 

In order to diagonalize Hq^ define the bare states (the arrows represent, from left to right, 
the z components of S*!, a and 5*2) : 

|i >= I TTi>, |2 >= I iTt>, |3 >= I TiT>, Sl^ = 1/2, (41) 

|4 >= I tii>, |5 >= I iiT>, |6 >= I iTi>, SI, = -1/2, (42) 

\up >= I TTT>, S^,, = 3/2. (43) 

Observe that S^^, = J2i S^ + a^ is conserved and thus each sector can be dealt with separately. 
Diagonahzing Hq in the sectors S^^, — ±1/2: 

det(< i\Ho\j > -e5i^) = (44) 

we get the eigenenergies Si which correspond to the eigenvectors: 

|i+>=X:a.,|j>,5^,, = l/2, (45) 

\i->=j2a,,\j>,Sl, = -l/2. (46) 

Both states are assumed to be normalized. The above two sectors have the lowest energy 
(compared to S^^, — ±3/2) and are degenerate. For definiteness we choose, for example 
S,^, — 1/2. It will be proven below, both analytically and numerically that the ground state 
of H stays in this sector for any L even after Hint is taken into account. 

In order to develop a perturbation theory in Hint it is convenient to introduce Green's 
functions. Since there are three states of three particles, we have to consider a 3 x 3 matrix 
three-particle Green's function Gij,i,j = 1,2,3. The poles of Gij are the matrix elements 
< i|i?o|j >■ In order to define Gij let us introduce a fermionic operator representation for 
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the three spins, analogously to (5). The fermions representing Si, S2 and a are and 
G^, respectively, /i =t, |. The Hamiltonian now should be thought of as expressed in terms 
of these operators. Then, for example Gii is defined as: 

Gii{t) = -i < r(*TW6»TW*iWMQT*t) >' (47) 
and similarly for all the others. The unperturbed Green's functions are: 

GT^\uj) = u- < i\Ho\j > +i6, z,j = 1,2,3. (48) 

To lowest order of perturbation theory with respect to Hint, the self-energies, corresponding 
to Gij are calculated as: 

= 1^1^ — - — + 



1=1 k 



£ — El — 



< i + \Hint\up, a >< up, a\Hint\j+ > ^^g^ 

]j £ ^up ^k 

where £„p — L/2 is the energy, corresponding to the state \up > (43). 

Then the equation for the effective energy level £* (compare with (30)) is: 

det{{£i - £*)6ij + ^ij{£*)) = 0. (50) 

The correct three-particle S^^^ — 1/2 eigenstate is given by: 

3 

|123 >= > . (51) 

i=l 

Here are the normalized eigenvectors of the matrix in (50). The three-particle - spin- wave 
ground state to lowest order is, similarly to (31): 

|G >= ^1 123 > +incoherent part. (52) 

The normalization factor Z is defined as: 

^=1+EW^^- (53) 
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The incoherent part is a mixture of the intermediate states > and \up,a > with 

appropriate weights, which are too lengthy to write down exphcitly. 

Finally, the average of any combination of spin operators O in this ground state can be 
easily calculated: 

< G\d\G >=< 123|6|123 >Z + Y^ \<up,a\H,nt\l23 > p ^ ^p\o\up > + 
3 < i23|iJ,„t|m-,/5></-,/3|if,„t|123> 

L L 7-— — Tvr^ — 7^ < - 1^1^- > • (^4) 

l,m=l k [e - Em - ek)[e - El - E]^) 

Using (54) we have calculated the spin-spin correlation functions for and a, as well as 
the appropriate magnetizations. The following notations are adopted: 

C((7, z) =< a.S, >, i = 1, 2, C(l, 2) S1.S2 >, (55) 

M(cr) a' >, M(i) =< 5f >, z = 1, 2 (56) 

We have also performed numerical simulations, based on exact diagonalization on small 
clusters with typical size N = 18 + 1. The numerical procedure is similar to the one outlined 
in Sec.IIIB. All our results are summarized in Figures 4 and 5. We have checked with both 
methods that the ground state is in the sector S^^^ = 1/2 for all L and thus no level crossing 
occurs. 

For ferromagnetic value of the coupling, L < 0, all correlations become ferromagnetic for 
sufficiently large |L|(Fig.4.). However, this is not accompanied by a change in the ground 
state, since all the three spins have a non-zero magnetization (Fig. 5.). This was already 
pointed out by Oitmaa, Betts, and Aydini. The situation is rather similar to the case of a 
FMB in an AFM background, studied in the previous section. 

For antiferromagnetic coupling, L > 0, there is a local ground state phase transition at 
L ^ 2.3, when the spins 5*2 and a change the direction of their magnetization. Naturally, the 
magnetization of Si increases at this point, since it is no longer frustrated. This local change 
in the ground state was observed numerically by Oitmaa, Betts, and Aydini. However, 
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they found a discontinuity in all correlation functions as well as the magnetizations at the 
transition point, which was attributed to level crossing. We, on the other hand, observe a 
continuous transition. It is quite possible that the discontinuity is due to a finite size effect. 

V. SUMMARY AND CONCLUSIONS. 

In summary, we have developed an analytic technique which allowed as to study the effect 
of locally frustrating perturbations (static spin defects) in 2D quantum antiferromagnets for 
any strength of the coupling. Such defects are of relevance to the physical system of oxygen 
holes in the CuO planes of the high-Tc cuprates at low doping. All analytic results are 
in good agreement with numerical simulations, based on exact diagonalization on small 
clusters. 

For an isolated ferromagnetic bond we found that the local magnetization does not vanish 
even for strong frustration. Our result suggests, that the previous calculations using lswaII 
lead to qualitatively wrong behavior. The region of validity of LSWA is thus restricted to 
small frustration only. 

We also studied frustration due to a quantum impurity spin coupled symmetrically to 
the two sub-lattices. It was found that for a ferromagnetic sign of the coupling there is 
no change in the ground state, similarly to the case of a FMB. For an antiferromagnetic 
interaction between the impurity and its neighbors we report a local change in the ground 
state. 

The technique presented in this paper is not limited to the study of the above problems. 
It can be applied to any spin systems where LRO is suppressed locally due to additional 
interactions and represents a nice alternative to LSWA for strong interactions. 
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FIG. 1. Diagrams for the self-energy S^: a) One-loop contribution. b),c) Examples of two-loop 
corrections. Solid lines represent fermionic Green's functions, while dashed lines are spin-wave 
propagators. 

FIG. 2. Local magnetization M (solid line) and the correlation function C(l,2) (long dashed 
line), calculated by using the Green's function method, as a function of the FMB strength K. The 
diamonds and squares are the corresponding exact diagonalization results. The short dashed line 
is the magnetization, calculated in LSWA. 

FIG. 3. The transverse (solid line) and the longitudinal (dashed line) parts of the spin-spin 
correlation function across the FMB. The triangles and the circles are the corresponding exact 
diagonalization results. 

FIG. 4. Correlations between the impurity spin and its neighbors, as defined in Eq.(55). The 
lines and the symbols represent, respectively, the analytical and the exact diagonalization results. 

FIG. 5. The magnetization of the impurity (solid line) and the neighboring spins (dashed lines) 
as a function of the interaction L. The symbols represent the exact diagonalization results. 
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